Interacting Arrays of Steps and Lines in Random Media 
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The phase diagram of two interacting planar arrays of directed lines in random media is obtained 
by a renormalization group analysis. The results are discussed in the contexts of the roughening 
of reconstructed crystal surfaces, and the pinning of flux line arrays in layered superconductors. 
Among the findings are a glassy flat phase with disordered domain structures, a novel second-order 
phase transition with continuously varying critical exponents, and the generic disappearance of the 
glassy "super-rough" phases found previously for a single array. 
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The statistical mechanics of planar arrays of directed 
lines is of interest to various physical systems. For ex- 
ample, the steps formed on a vicinal crystal surface can 
be modeled by an array of directed lines confined in a 
plane @. The same model describes Josephson vortex 
lines in a planar Josephson junction subject to a parallel 
magnetic field . An important issue concerns the be- 
havior of such arrays of lines in the presence of quenched 
disorders, ft is known that the lines are pinned by point 
disorders and become "glassy" at low temperature [Q. 
However, different analytic and numerical || stud- 
ies have yielded conflicting results regarding the nature 
of the glass phase. A subject of debate is the possib le 
breaking of replica-symmetry in the glass-phase 



possibK 




FIG 1: Two kinds of (3x1) microfacets on a (2x1) re- 
constructed crystal surface. The background (2x1) facets 
can be on four possible sublattices (marked "a"-"d"). 
Each (3x1) facet shifts the phase by one sublattice. 

In this article, we study the effect of point disorders 
on two interacting species of lines in planar arrays. This 
problem arises in the study of anisotropically (2x1) re- 
constructed gold (110) surfaces where two kinds of 
(3x1) microfacets can be treated as two species of inter- 
acting lines [Fig. 1]. Previous studies of the pure sys- 
tem have revealed a rich phase diagram with a variety of 
possible phases as a function of the interaction parame- 
ters [ ]l2| , |l3"| | . The inclusion of point disorders, say crys- 
talline defects originating from a disordered underlying 
substrate, induces deformations in the trajectories of the 
microfacets. Similar issues arise in two layers of magnet- 
ically interacting Josephson vortex lines. Performing a 



renormalization group (RG) analysis in replica space, we 
are able to obtain a complete picture of the RG-flow. The 
result is applied to discuss the rich phase diagrams of the 
anisotropically (2x1) reconstructed crystal surfaces. We 
also discuss the structure of the glass phases obtained for 
two planar vortex arrays, and comment on the relevance 
of our results to the issue of replica-symmetry breaking 
in a single vortex array. 

A single species of directed lines confined in a plane 
containing quenched randomness can be described by the 
continuum Hamiltonian 



y(r)p(0(r),r) (1) 



on length scales exceeding the line spacing I. The first 
part of (P gives the elastic energy of the line array in 
terms of a displacement-like scalar field </>(r) (a displace- 
ment by I corresponds to a shift of 2tt in <$), charac- 
terized by an (isotropized) elastic constant K [M. The 
second term describes density variations p(4>, r) induced 
by a random potential V(r). The density field has the 
form p(<f>(r), r) w p (l- d x cj>/ (27rp ) + 2 cos [2-Kp x -<£]), 
where p = l/l is the average line density, and r = (x, z) 
with z along the line direction. The random potential 
is taken to have zero mean with short-range correlations 
V(r)V(r') — gS 2 (r — r') of (bare) strength g. 

An interaction between two such species of lines in the 
form of / ri)1 . 2 Vint(ri-r 2 ) / o(0i,ri)p(^2,r 2 ) with a short- 
ranged potential Vi nt p5[ leads to 



(3H 
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2fip Q cos ( 
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with p = J r Vi n t(r) and = p/8n 2 Q. We assume the 
disorder potential Vt acting on species i to be statisti- 
cally identical, with the cross-correlations Vi(r)V2(r') = 
g^5 2 {v — r') to be specified shortly. The full Hamil- 
tonian of our system, H[<f>i,<j>2] = Yli=i ^ D > ^ + 
Wint [</>].) 4>2\, can then be written in a succinct form (af- 
ter neglecting irrelevant terms |14|), 
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r i,j i 

-J2W i (<t>i,r)+2fxp? cos((f> 1 -<j> 2 )\, (3) 

i ^ 

with effective random potentials Wi and whose corre- 
lators are Wj(0, r)W 3 (0', r') = 2g i] p\ cos (0—0') <5(r-r'), 
and w i (r)w :) (r') = A i:) 5(r — r'), with the latter having 
the bare values Ay = gij/(8w 2 ). The parameters of the 
theory are — K,gij — g,A.ij = A for i = j and 
Kij = K^, gij = g^, Ay = A M for i ^ j. The appearance 
of the cosine couplings reflects the discrete translational 
symmetry of the line arrays and has been extensively dis- 
cussed in the literature. 

A physical observable of interest for the reconstructed 
crystal surface is the interface height profile h(r). It 
is given by h(r) = [(f>i(r) — 02 (r)] /2ir since lines from 
species 1 and 2 represent upward and downward steps 
respectively, and it is the difference of the two that de- 
termines the height profile (see also Ref. [|l3|). Note also 
that the (2x1) reconstructed crystal has a choice of four 
possible sublattices on the surface (see Fig. 1). Each step 
is thus also a domain wall separating domains of the sub- 
lattices. An order parameter capturing the ordering of 
the domains is given by R(r) = exp[i(0i(r) + 02(r))/4]. 
To find the large scale behaviors of these observables, 
we apply the Replica method to average over disorders. 
A replica-symmetric RG analysis fl6| for the replicated 
system yields the following recursion relations (to bilin- 
ear order) upon a change of scale by a factor e l : 

dg/dl = (k- r)g - g 2 - g^ , 
dg^/dl = (k-t- 8)g ll - gg^ - g[i , 
dp/dl = (2k — 6) (j, — (t + n)g^ — g\i , 
dn/dl = fx(fi - 2 ffM )/2, dS/dl = (g 2 - fi)/2 , (4) 

where k = 1 — [Air(K — K^)]^ 1 is a reduced elasticity 
parameter, 5 — 87r(A — A^), and a non-universal numer- 
ical factor is absorbed into g, \i and g^. The RG flow 
is controlled by the reduced temperature, r = [A-n(K + 
K^)]^ 1 — 1, which is not renormalized due to a statis- 
tical tilt symmetry. [There exists a sixth RG equation 
d5/dl = (g 2 + gf t )/2, for the parameter 5 = 8tt(A + A M ). 
It however does not feedback into (|J).] 

Before delving into the structure of the RG flow, we 
first mention two limiting sub-problems which have been 
studied previously. In the limit fi, g^ = 0, the elastic- 
ity parameter n is also not renormalized, and the RG 
equation has the same structure as that obtained for a 
single species of lines Q , with an effective temperature 
T c ff = t — k. As shown in Fig. 2 (inlet (a)), The disorder 
(g) is irrelevant at high temperatures (r > k), yielding 
the usual logarithmic roughness for 2D surfaces, accom- 
panied by a quasi-long-ranged domain order, with power 
law decays in (R(r)R(r')) . We refer to this as the decou- 
pled line (DL) phase. At r = k, the marginal irrelevance 



of g yields a marginally-coupled line phase (ML), which 
again has logarithmic roughness and quasi-long-ranged 
domain ordering. At low temperatures (r < k), the dis- 
order is relevant. The resulting glass phases are described 
by the line of fixed points g* = n — r which are pertur- 
batively accessible for \k\ , |r| <C 1. There, one finds the 
interface to be super rough Q, with (h 2 ) ~ log 2 L (due 
to the divergence of 8) on large scales L. The glassiness 
is also reflected by a disordered (short-ranged) domain 
order, i.e., a faster than power-law decay in (R(r)R(r')) , 
due to the divergence of 5. Complete decoupling for the 
interaction (||) implies that = also, or n + t = 0. 
Thus g* = 2\t\ (and fi* = 0, k* = |r|) is the physi- 
cal fixed point describing a decoupled glass (DG) phase. 
Another well known limit of our problem is that of van- 
ishing disorder {g,g^ = 0), where a Kosterlitz-Thouless 
(KT) transition occurs independent of r (see Fig. 2 in- 
let (b)). For large coupling fi, the two species become 
locked together, forming an elastically-coupled line (EL) 
phase with fi* , k* — > O(l) (and g* — 0). Since the up 
and down steps are now paired, the interface is flat, with 
again quasi-long-ranged domain order. 




FIG 2: RG-flow for identical disorder (g^ = g, A = 
A M ). Inlets (a) and (b) describe the flow for the sub-case 
with no inter-species coupling (p = 0) and no disorders 
(g = 0), respectively. 

Our task is to obtain the flow diagram of the full sys- 
tem with nonzero (but initially small) disorder g and in- 
teraction p. We shall focus on two choices of disorder 
which are of particular interest: (i) identical disorder for 
the two species (g^ = g, S = 0), and (ii) completely 
uncorrelated disorder (g^ — 0) . The case of identical dis- 
order is the generic situation for steps on reconstructed 
crystal surfaces with disorder. It can also be constructed 
for the two layers of Josephson vortex lines. The most 
striking feature of the RG flow in the low temperature 
(r < 0) regime [Fig. 2] is the strong instability of the 
DG fixed point with respect to inter-species interaction 
/i 7^ 0. An attractive interaction (fi < 0) favors the two 
species to lock into the same configuration, i.e. 0i = 02- 
Once locked, the system acts effectively as a single species 
with a doubled elastic constant K, or equivalently a lower 
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effective temperature T c ff < r, so that the effective single 
species problem is in the glass phase. Correspondingly, 
(||) yields an RG flow away from the unstable DG fixed 
point to a sink with strong inter-species coupling and 
strong disorder [g*,K*,-^* — > 0(1)], for all r < 0. We 
refer to this as the elastically-coupled glass (EG) phase. 
Because opposite steps are tightly bound, the EG phase 
is flat. However, the domain structure is disordered (as 
in DG) due to the pinning of the bound pairs of steps (or 
domain walls). For r > 0, we find the RG-flow to be dom- 
inated by a KT-like transition in the parameters (k,/i), 
similar to that of the disorder-free sub-problem, with the 
strongly-coupled sink being the low temperature EG. 

A repulsive interaction (/i > 0) competes with fluc- 
tuations in the random potential, which still attempts 
to lock the two species into the same configuration. At 
low temperatures (t < 0), this competition leads to two 
RG sinks separated by a second-order phase transition: 
If the repulsive interaction dominates, the two species 
avoid each other by locking into a configuration with 
4>i = (p2 + 7r, i.e., with one species displaced by half a 
line spacing with respect to the other. Such a configu- 
ration can be interpreted again as a single species, but 
now with a doubled line density. This leads to a higher 
effective temperature T e s > 0, such that the effective 
single species system is not glassy, with g* — 0. Corre- 
sponding to this scenario, we find for weak bare disorder 
a RG-flow away from DG towards the fixed point EL, 
and the RG trajectories approach their pendants in the 
disorder-free sub-problem. If disorder dominates, how- 
ever, we obtain a RG flow from DG towards the fixed 
point ML, since the disorder weakens the inter-species 
coupling [i and k, while the coupling fj, in turn weak- 
ens the disorder g. The phase transition separating the 
EL-phase and the ML-phase is governed by an unstable 
fixed point (T) at (K*,g*,n*) = (-1, 2, 4) \r\ /7, which 
is the attractor of the plane of separatrix g ~ /i + 2k. 
The phase transition is second order for r < 0, with an 
algebraically diverging correlation length (characterized 
by an exponent v — 2 |r| /7) upon crossing the separa- 
trix. It is a remarkable feature of this system that while 
both RG sinks are (at least marginally) disorder- free, the 
unstable fixed point governing the transition is disorder- 
dominated, as manifested by g* > and a disordered 
domain structure. 

The above analysis of the RG flow can be straight- 
forwardly turned into a phase diagram. We define an 
inter-species interaction energy U = fiT, and present 
the phase diagram in the (U,T) space for a fixed dis- 
order (g) [see Figs. 3]. The super- rough DG phases 
exist at U = below a critical temperature T c given 
by K c = 1/ (47r) . These phases are marked by the 
thick wavy line in Fig. 3(a). With attractive interaction 
(U < 0), there exists at low temperature a flat, disor- 
dered phase (EG) where opposite steps are bound and 
pinned by disorders, with a KT-like transition (dashed 



line in Fig. 3(a)) to a logarithmically-rough, pure phase 
(DL) where steps unbind. A rather different scenario is 
obtained for repulsive interactions (U > 0): The separa- 
trix g = [i + 2k separates the flat, pure phase (EL) at low 
temperature and large repulsion from the two logarithmic 
phases. At very high temperatures (T ^> T c ), the sys- 
tem is in the pure decoupled phase (DL), with the usual 
logarithmic interface roughness and quasi-long-range do- 
main order. Upon lowering the temperature beyond the 
line t = (thin solid line), the system settles into the 
ML phase for weak repulsive interaction (compared to 
the strength of disorders). The ML phase is in practice 
quite similar to the high temperature phase, except for 
differences in the coefficient of the logarithmic interface 
roughness and the exponent of the algebraic decay of do- 
main order. Further lowering the temperature beyond 
the separatrix (the thick solid line), the system makes a 
second-order transition from the ML to the EL phase. 
Note that because the critical properties there are con- 
trolled by the fixed point T which depends on r, the 
critical exponents governing this transition actually vary 
continuously along the thick solid line. The second-order 
transition terminates at a point where the separatrix in- 
tersects the line r = (the open circle in Fig. 3(a)). The 
transition between DL and EL at higher temperatures 
still occur along the separatrix (dashed line) and is of 
the KT-type. 




FIG 3: Phase diagram for two species of lines in (a) 
identical and (b) uncorrelated random media. 

A curious result of the above analysis is the generic 
disappearance of the glass order (as reflected by super- 
roughening and/or disordered domain structures) upon 
the introduction of repulsive interactions. Such glass or- 
der was found previously for a single line array in random 
media using the replica-symmetric RG method Jl],[|,|]], 
but was absent in more recent studies using a variational 
method with replica-symmetry breaking (RSB) [pl[|. 
The latter, which finds instead the ML phase, appears 
naively consistent with our findings. On the other hand, 
our result can be used to question the internal consis- 
tency of the RSB scheme: As described in Ref. fnj, a 
physical way of probing the existence of RSB is to take 
two physical replicas of a system (in identical random 
potential), and monitor the response to a small repul- 
sion between the replicas. If there is a degeneracy of 
low free energy states (which the RSB scheme attempts 
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to describe), then an infinitesimal repulsion between the 
replicas will force the two to occupy different states which 
have similar glassy properties and little overlap. The sys- 
tem we have analyzed so far can be interpreted as two 
physical replicas in the same random potential. In our 
case, we note that the repulsion has so strong an effect 
that it yields completely different glassy properties, i.e., 
from DG to ML. This indicates that the glass order for 
the single line array is extremely fragile |l^,|ll|, making 
it quite different from the usual scenario expected of sta- 
ble (zero-temperature) glass phases. Thus from the view 
point of the replica-symmetric RG analysis, the absence 
of glass order from the solution of the variational treat- 
ment is not surprising, as it may be the result of subtle 
interactions introduced by the RSB scheme itself. 

We continue with a short discussion of the case where 
the disorder potentials acting on the two species are un- 
corrected, i.e. with the bare g^ = 0. This is the generic 
case for two layers of Josephson vortex lines in disor- 
dered planar Josephson junctions. Uncorrelated disorder 
tends to decouple the two species of lines and competes 
with the locking effect of the inter-species interaction. 
(Here attractive and repulsive interactions are qualita- 
tively similar, up to a relative 7r-phase shift between the 
two species.) The RG analysis is more complicated than 
before, because one has to consider the full set of equa- 
tions in (^). [Note that g^ is generated by \x and g.] The 
basic features of the phase diagram can be obtained by 
observing that the RG-flow is dominated by two KT tran- 
sitions which can be found in two sub- problems of (Q): 

(i) The KT transition of the disorder-free sub-problem. 

(ii) The KT transition which occurs in the space (g^, S), 
when g = fi = k + t = 0. The latter has a critical separa- 
trix g^/^/2 = 8 + 2r which isolates two regions of flow to 
a sink with g^ — and 5 > 0, and another sink where g^ 
grows. The first sink is consistent with a decoupled glass 
(DG) phase as the eigenvalue for the flow of /i becomes 
negative there, while the second sink is consistent with 
the elastically-coupled glass (EG) phase where 1^1 — > g. 
The form of this separatrix for the regime of physical 
initial conditions (e.g., g,pL^> g^) is not known analyti- 
cally. However, we determined it numerically to be well 
described by the form 



\i ~ eg, 



(5) 



with the numerical constant c w 0.5. This condition, 
combined with with sub-problem (i) leads to the phase 
diagram depicted in Fig. 3(b). It is interesting to note 
that the condition (^) separating the decoupled and cou- 
pled glass phases is very similar, including the numerical 
value of c, to the result obtained recently for a problem 
involving many layers of vortex lines |1S| ] . 

In conclusion, we have presented a detailed RG anal- 
ysis for a model of two interacting planar line arrays in 
random media. Among the findings are a novel second- 
order phase transition with continuously varying critical 



exponents, a disordered flat phase, and the replacement 
of the super-rough glass phase by a marginally-coupled 
phase. These findings along with the structure of the 
proposed phase diagrams should be accessible by exper- 
imental or numerical investigations. 
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Note added: After the completion of this work, we re- 
ceived a preprint by Carpentier et al Jl9[| who obtained 
the same RG equations for the two species of interacting 
lines as our Eq. |J. The analysis for the case of uncor- 
related disorders (<? M = 0) described in Ref. |L9| agrees 
qualitatively with our finding. The case of identical dis- 
orders^ = g and 6 — 0) was however not discussed in 
Ref. [fll. 
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